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Motivated by the recent experiment of Wang et al. [Nature Physics 6, 389 (2010)], who observed 
a highly unusual transport behavior of ferromagnetic Cobalt nanowires proximity-coupled to su- 
perconducting electrodes, we study proximity effect and temperature-dependent transport in such 
a mesoscopic hybrid structure. It is assumed that the asymmetry in the tunneling barrier gives 
rise to the Rashba spin-orbit-coupling in the barrier that enables induced p-wave superconductivity 
in the ferromagnet to exist. We first develop a microscopic theory of Andreev scattering at the 
spin-orbit-coupled interface, derive a set of self-consistent boundary conditions, and find an expres- 
sion for the p-wave minigap in terms of the microscopic parameters of the contact. Second, we 
study temperature-dependence of the resistance near the superconducting transition and find that 
it should generally feature a fluctuation-induced peak. The upturn in resistance is related to the 
suppression of the single-particle density of states due to the formation of fluctuating pairs, whose 
tunneling is suppressed. In conclusion, we discuss this and related setups involving ferromagnetic 
nanowires in the context of one-dimensional topological superconductors. It is argued that to realize 
unpaired end Majorana modes, one does not necessarily need a half-metallic state, but a partial spin 
polarization may suffice. Finally, we propose yet another related class of material systems - ferro- 
magnetic semiconductor wires coupled to ferromagnetic superconductors - where direct realization 
of Kitaev-Majorana model should be especially straightforward. 

PACS numbers: 74.40.+k, 74.81. Bd, 64.60.Ak 



I. INTRODUCTION 

Coexistence of ferromagnetism and superconductiv- 
ity is a rare phenomenon, yet there exists unambigu- 
ous evidence for its occurrence. Superconductivity 
has been thoroughly investigated in the uranium-based 
itinerant ferromagnets^"^ and there have been reports 
about possible co-existence of these phases in the d- 
electron compound ZrZnJ^ and the copper oxide com- 
pound RuSr 2 GdCu 2 8 _ 5 (Ru-1212^1. The neutron 
diffraction measurements in ferromagnetic superconduc- 
tors ErRh 4 Bjin] HoMogSj 111 and HoMo 6 Se^2l indicate 
an inhomogeneous magnetic order coexisting with su- 
perconductivity. More recently, coexistence of these 
seemingly exclusive states was also seen in Pb/PbO 
core/shell nanoparticles^, and in two-dimensional inter- 
faces betwe en perovskite band insulators LaA103 and 
SrTiOjJSE!] 

The interplay of superconductivity and ferromag- 
netism can also be studied in the context of supercon- 
ducting proximity effect, where a superconductor in con- 
tact with a normal metal induces superconducting cor- 
relations in the latter. Due to incompatible spin or- 
der conventional superconducting correlations are known 
to penetrate negligibly inside a ferromagnet However, 
a number of recent studies have demonstrated an un- 
expectedly long-ranged proximity effect in mesoscopic 
superconductor-ferromagnet hybrid structures EHH1 A 
subsequent theoretical work showed that local inhomo- 



geneity in the magnetization near the interface can in- 
duce triplet-pairing inside the ferromagnet, and that 
these correlations can account for the observed long- 
ranged proximity effect This beautiful theory predicts 
an exotic odd-frequency (even-momentum) symmetry for 
the triplet component of the condensate, which was orig- 
inally suggested by Berenzinskii as a possible phase in 
superfluid 3 He Theoretical work on the Josephson cou- 
pling between two conventional superconductors sepa- 
rated by a half-metallic ferromagnet has also investigated 
the mechanism behind singlet-to-trip let co nversion near 
a half-metal-superconductor interface! 25 * 26 ! 

Seemingly unrelated to the developments in supercon- 
ducting proximity effect at the time, Kitaev showed in 
his pioneering work that Majorana fermion excitations 
can be localized at the ends of a spinless p x + ip y su- 
perconducting quantum wireP^l Conceptually, the Kitaev 
model of a topological superconductor is that of a fully- 
polarized ferromagnetic superconductor. However, a di- 
rect connection between the model and the existing ferro- 
magnetic itinerant and hybrid superconducting systems 
has been made only very recently. In some of the re- 
cent proposals, superconducting proximity effect plays a 
key role in realizing topological superconductors which 
support Majorana fermions on their boundaries or in a 
vortex corePSES 

On the experimental front, Ref. [21] has already shown 
that a ferromagnetic cobalt (Co) quantum wire can be 
made superconducting by placing it in contact with a 



conventional superconductor. Proximity-induced super- 
conductivity was observed in the wire over a distance 
of a few hundred nanometers. The pairing symmetry 
inside the wire is likely to be a spin-triplet state since 
singlet correlations should not survive inside the fer- 
romagnet over a distance exceeding a few nanometers. 
While the odd-frequency superconductivity indeed re- 
mains a viable explanation of the long-range proximity 
effect, an alternative scenario in clean wires involves p- 
wave superconductivity induced by a spin-orbit coupled 
interface.^ The latter would imply that an experiment 
of the kind in Ref. [5T]is very close to realizing the topo- 
logical superconductor^, and may host the sought-after 
Majorana fermions at the two ends of the wire or inside 
the wire, wherever the parity of the number of occupied 
sub-bands changes. 

With its broader physical relevance aside, the experi- 
mental work of Ref. [21] has made very interesting obser- 
vations in the context of superconducting proximity effect 
in ferromagnetic systems. The work systematically stud- 
ies the resistance of ferromagnetic single-crystal cobalt 
nanowires sandwiched between two superconducting elec- 
trodes. They observe proximity-induced superconductiv- 
ity in the nanowire over a distance of order 500nm, which 
is orders of magnitude longer than the coherence length 
expected for conventional superconducting correlations 
inside a ferromagnet. For some wires, transition to super- 
conductivity is preempted by a large and sharp resistance 
peak near the transition temperature of the electrodes. 
The resistance peak disappears when the Co wire is re- 
placed by a gold wire, indicating an intimate connection 
between it and the ferromagnetism of the wire. 

Motivated by these observations, we first theoreti- 
cally study the physics of an s-wave-superconductor- 
ferromagnet tunnel junction. We assume that the asym- 
metry in the tunneling barrier gives rise to the Rashba 
spin-orbit coupling inside the tunnel barrier and en- 
ables p-wave superconductivity to be induced inside the 
ferromagnet We derive a complete set of relevant 
boundary conditions, and an expression for the p-wave 
minigap is found in terms of the microscopic parame- 
ters of the contact. In the second part of the work, we 
investigate in detail how the anomalous resistance peak 
can be explained based on the theory of supe rconducting 
fluctuation corrections to conductivity!^ 1 ^! We discuss 
a possible faithful model for the cobalt nanowire in the 
vicinity of the deposited W electrodes, develop a micro- 
scopic theory for the superconducting fluctuations in the 
W electrodes, and study how these fluctuations influence 
the resistance of the cobalt wire. Finally, we discuss in 
more detail the possible relevance of the experiment to 
topological superconductivity and propose a related hy- 
brid system, which provides in our opinion the simplest 
physical realization of topological superconductivity and 
localized Majorana modes. 

We emphasize that the Co wire throughout the work 
is treated mostly in the ballistic/ballistic-to-diffusive 
crossover limit. As we outline in later parts of the paper, 



the treatment of the wire in this limit is reasonably jus- 
tified based on experimental datepH and past studies on 
Co band structure^lED although the information avail- 
able to us is insufficient to conclusively determine the 
nature of proximity-induced superconductivity in the hy- 
brid system. As we have mentioned above, long-ranged 
odd-frequency p-wave proxi mity e ffect, extensively stud- 
ied in diffusive ferromagnet sr 2 ^ 23 !, remains another real- 
istic scenario. Here, we are exploring the possibility of 
a long-ranged spatially- odd p-wave proximity effect in- 
duced via intcrfacial spin-orbit coupling, which remains 
a viable physical scenario when the ferromagnet is in the 
clean limit. Furthermore, as argued below, this latter 
mechanism may survive even in disordered systems due 
to mesoscopic fluctuation effects. 

The paper is organized as follows. Sec. |Tl| summarizes 
the main findings of the experiment in Ref. I2R with a 
particular emphasis on the critical resistance peak. The 
section also lays out in detail an estimate for the mean- 
free path of the Co wire and a justification for consid- 
ering the clean limit. A microscopic theory for Andreev 
scattering at a ferromagnet-superconductor interface is 
presented in Sec. |III| In Sec. |IV| the superconducting 
fluctuation theory is applied to the W electrodes. Fol- 
lowing a qualitative discussion on how superconducting 
fluctuations in the W can lead to resistance peaks in Co 
transport, a detailed microscopic theory for the correc- 
tions to resistance from these fluctuations is developed in 
Sec. |IV A| A discussion of topological superconductivity 
is presented in Sec. |V| 



II. SUMMARY OF THE PENN STATE 
EXPERIMENT 



Ref. [3T] reports observations of long-ranged proximity 
effect in single-crystal ferromagnetic Co nanowires. Re- 
sistance of the wire is studied using a four-probe setup 
with all four electrodes made from superconducting tung- 
sten (W) [see Fig. TJa)] or with a combination of super- 
conducting W and non-superconducting platinum (Pt) 
electrodes [see Fig. [TJb)] . While the outer electrodes are 
used to pass current i w through the wire, the inner elec- 
trodes are used to measure the potential difference across 
a length L of the wire. With superconducting voltage 
electrodes [Fig. [l|a)] Co wires become fully supercon- 
ducting below the electrode transition temperature for 
L = 600nm, and show substantial resistance drop for 
L = 1.5 //m. However, such resistance drop was not ob- 
served when the superconducting electrodes are replaced 
by non-superconducting Pt. However, when an addi- 
tional electrically isolated W strip is deposited between 
the Pt electrodes (see Fig. [l|c)), proximity-induced su- 
perconducting properties are restored in the Co wire. 
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FIG. 1: A cartoon representation of the various experimen- 
tal setups used in Ref. 1211 Three main setups are considered: 
(a) W current and voltage electrodes; (b) W current elec- 
trodes and Pt voltage electrodes; and (c) same as (b) with an 
additional electrically isolated W strip between the Pt elec- 
trodes. Length of the wire L is determined by the distance 
between the inner edges of the voltage electrodes. The width 
(or diameter) of the Co wire is denoted by w. 



A. Critical resistance peak 

The most interesting observation in the experiment is 
that the resistance shows a sharp peak as a function of 
temperature as it approaches the superconducting tran- 
sition temperature of the W electrodes. A cartoon pic- 
ture of a typical resistance versus temperature curve is 
shown in Fig. [2] to illustrate this result. The transition 
temperature of the electrodes is estimated to be between 
T c = 4.4-5K. This resistance peak is very large and con- 
stitutes 25-100% of the normal state resistance depend- 
ing on the wire width. Intriguingly, the peak disappears 
when the Co wire is replaced by a paramagnetic gold 
wire, and thus seems to be a consequence of the ferro- 
magnetism of the Co wire. The large resistance peak 
also disappears when the W electrodes are replaced by 
Pt electrodes [see Fig. [ljb)], but is restored once an elec- 
trically isolated W strip is deposited [see Fig. [TJc)]. The 
large resistance upturn thus seems to require only that 
a superconducting strip is in contact with the nanowire, 
and does not require it to be either a current or a voltage 
electrode. 
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FIG. 2: A cartoon plot of the experimentally observed Co 
wire resistance as a function of temperature. A critical resis- 
tance peak is observed near the transition temperature T c of 
the W electrodes. A precipitous drop in the resistance is then 
observed below T c . 



B. Mean-free path 

When the extra W strip is deposited on the wire [see 
Fig. [ijc)] data show that the wire's normal state resis- 
tance increases by approximately 55f2. This amounts to 
a nearly 50% increase from its normal state resistance 
of «126f2 prior to the deposition [see Fig. [ljb)]. This 
suggests that the deposition process could have strongly 
modified the property of the wire in the vicinity of the 
contact. Returning to the setup in Fig. [lja) the exper- 
imental data show total normal state wire resistance of 
approximately 145S1. The resistance of the Co wire aris- 
ing from regions unaffected by the W strips can then be 
approximated as 145f2 — 2 • 55J7 = 3517. Using the quoted 
distance between the voltage electrodes (i.e. L = 1.5/xm) 
and the wire width (i.e. w = 36nm) the corresponding 
resistivity is p ~ 23nSl-m. An estimate for Co's Fermi 
velocity, based on critical current oscillations observed in 
a Nb/Co/Nb Josephson junction, gives Dp ps 280km/sP^ 
From this the density for electrons with a single spin pro- 
jection can be estimated as 



3„,3 



, vrctv 



6n 2 h 3 



2.4 x 10 26 m 



(1) 



where x = m* /m e is the Co effective mass ratio with m* 
being the effective mass of the Co electrons and m e the 
mass of the electron. The Drude formula then gives a 
mean-free path estimate of 



xm e VF 
ne 2 p 



1.8/xm 



(2) 



Cores of the Co nanowires studied in Ref. [21] show a 
hexagonal close-packed zone pattern and the wires have a 
[0001] growth direction. Presumably, the electron trans- 
port through the wire is then along the c-axis direction. 
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For our estimate of the Co effective mass x, we use results 
based on band structure calculations for Co cyclotron 
mass.^ In the relevant direction of transport, we have 
x « 3. Inserting this value into (2| the mean- free path 
is estimated to be t ss 200nm, which is the same order 
of magnitude as the observed coherence length in the Co 
wire. Within this scenario, regions of the Co wire unaf- 
fected by the W electrodes can be considered in the clean 
limit. 

We note, however, that the above estimate for the 
mean-free path relies crucially on our estimate for the 
resistance of the Co wire unaffected by the W electrodes. 
Furthermore, the estimate also depends sensitively on the 
values used for vp and x, which we do not know with 
certainty. Due to Co's complicated band structure, the 
values for Vp and x strongly depend on the direction of 
transport with respect to the crystallographic axes and 
on the Fermi surface which gives the dominant contribu- 
tion to transport. For instance, for vp ~ 10 6 m/s (quoted 
in Ref. I2T|) we obtain a much shorter mean-free path 
of £ f=a 16nm, implying that the wire is in the diffusive 
limit. Co having multiple subbands with varying x also 
implies that some subbands are less affected by disorder 
than others. If all bands participating in transport indeed 
have short mean- free paths (of order a few nanometers), 
the likely mechanism for pr oximit y effect is the spatially 
even, odd-frequency pairing. 22 2 - However, based on ex- 
perimental data in Ref. [2]] and on Co band structure and 
de-Haas-van Alphen studies relevant for electron trans- 
port in the c-axis direction^, we have shown above that 
there is reasonable justification to model the Co wire in 
the clean limit, and the spatially odd p-wave pairing is a 
viable mechanism for the observed proximity effect. Fur- 
thermore, mesoscopic fluctuation effects may further en- 
able remnants of p- wave superconductivity to persist well 
into the diffusive limit In this work, we pursue this 
scenario of a Co wire with proximity-induced p-wave su- 
perconductivity, and purport that it provides a natural 
explanation for the resistance peak. 



III. ANDREEV REFLECTION AT A FM-SC 
INTERFACE 

We consider a tunnel junction consisting of a ferro- 
magnetic normal metal at < z < d, a tunneling barrier 
at —a < z < 0, and a superconductor at z < —a (see 
Fig. [3]). At the moment, we assume translational sym- 
metry in the direction(s) parallel to the interfaces. We 
will work with a planar junction, but the results will also 
describe proximity effect in a ferromagnetic nanowire by 
simply "removing" one of the two transverse directions or 
by choosing a specific transverse sub-band. The deriva- 
tion below follows closely that of Refs. |49] and [50l but 
is generalized to include spin-orbit-coupling effects in the 
barrier. 
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FIG. 3: The proximity tunnel junction considered. The 
ferromagnetic normal metal, the tunnel barrier with spin- 
orbit coupling and the superconductor occupy < z < d, 
—a < z < and 2 < —a, respectively. 



A. Equations of motion 

The electron field operators in the ferromagnet ip^ F \ 
tunnel barrier and superconductor ?/j( s ) satisfy the 
following equations of motion. In the normal region (0 < 
z < d), 



H2d2 A x 



O)=0; (3) 



in the barrier (—a < z < 0), 



) 2 - h 2 d 2 z 

2m B 



5 a - a R [p x er Q/3 ]j 



x4%)=0; (4) 
and in the superconductor (z < —a), 

( £ -^ s)+ S)^ (z) " A ^ SHp(z)=o - (5) 

Here, p = (p x ,P v ) is the two-dimensional momentum 
parallel to the interface, <x = (a x , (J yi & z ) is the vector of 

Pauli matrices acting in spin space, £^ — p 2 /2mp — \ip 

and — p 2 /2ms — ar e the electronic spectra in 
the ferromagnet and superconductor, mp, ms and tub 
are the effective masses in the corresponding region and 
lip and us are the Fermi levels in the ferromagnet and 
superconductor, a and f3 here label the spin projections. 
We assume that the ferromagnet and superconductor are 
separated by a tunnel barrier of height Uq that, due to 
mirror asymmetry with respect to reflections relative to 
the z = —a/2 plane, contains Rashba spin-orbit cou- 
pling characterized by parameter a^. The ferromagnet 
is modeled as a Fermi gas in a Zeeman field h. Finally, 
the superconductor is assumed to be of s-wave type and, 
consequently, the mean-field order parameter is momen- 
tum independent and reads A Q ^ = A (ia y, a p)- 

Although the barrier is spin-orbit coupled, this cou- 
pling will not lead to any spin Hall effects at the bound- 



4 



aries z = and z = — a. Consequently, there is no accu- 
mulation of spin and the basic ballistic boundary condi- 
tions are identical to those in the theory of a conventional 
tunnel junction, namely: 



0, 



(6) 



$S(0) = O), ^(-«) = ^(-a), (7) 



-5 2 ^(0) = —9^(0), 



and 



m 5 



(8) 



(9) 



The full set of equations (|3|-([5]) is formulated for elec- 
tron field operators acting in Fock space. However, since 
these equations and boundary conditions (|6|-(|9| are all 
linear, the operator nature of the unknown functions is 
not germane and we may simply approach the problem 
as in single-particle quantum mechanics. The problem 
is then conceptually very simple and reduces to solving 
second-order differential equations albeit with non-trivial 
boundary conditions. This can be done exactly to a large 
degree. 



B. Deriving closed boundary conditions 
for the ferromagnet 

1. Tunnel barrier 

Let us look for a solution in the tunnel barrier in the 
following form, ipa : p(z) — Z(z)V a (p), where the trans- 
verse wave function Z(z) satisfies 

h 2 „, 



2m f 



U -e + e H (p)]Z = 0, 



(10) 



and the planar wave function V a (p) is a standard wave 
function of the Rashba problem, 



^(P) 



1 



±1 



(11) 



describing the helicity eigenstates with eigenvalues, 
£±(p) = 2ST^ ^ a fl|p|- The angle 7 P is defined via 
p = p (cos 7 p , sin 7 p ) . The general solution in the bar- 
rier then reads 



[C- 



■C-e-«-*]V-(p), (12) 



2mB [U - £ + e±(p)]. To guarantee real 



where q± — y 
q± we assume that the spin-orbit interaction scale is 
smaller than the barrier height. 



We now use ( 7| to determine the relations between the 
coefficients in (12) and the values of the wave function 



on the superconductor and ferromagnet boundaries. We 
find 



Civ 



(jrl e l ' q,a -S l ) 
,a) V J 



4 sinh(<7;a) 
where 1,1' — ± and 

^=#2(0)e-^±^<2(0); 
5± -.-./,<*)/ 



i^>(-a)e-^±^(-a) 



(13) 

(14) 
(15) 



2. Superconductor 

We now match the tunnel-barrier solution ( fl2| ) 
with the superconducting solution at the barrier- 
superconductor interface using boundary condition 
We then find 

<p(-«) = r p , a ^S(0) - flp,a^S(-0)- (16) 



The matrices in (16), which play an important role in the 



analysis of the Andreev scattering problem, are given by 



T P = 



«f 



k i8ne * 7p 
—ibne llv k 



Here, we have defined 



and 



8k = 



q+ + 


Q- 




2 sinh(q + a) 
<?+ 


2 sinh(g_ 

q- 


a) 


2 sinh(q + a) 


2 sinh(g_ 


a) 


q+ + 


q- 




2 tanh(g + a) 

q+ 


2 tanh(q_ 

q- 


-a) 



2 tanh(<7 + a) 2 tanh(g_o) 



(17) 
(18) 

(19) 
(20) 

(21) 
(22) 



with q± l being the penetration length of a particle with 
a positive/negative chirality inside the barrier, and a is 
the barrier width. As the width of the barrier grows, the 
"tunneling" boundary coefficients decay exponentially, 
lim a _ i . 00 n t — 0. As expected, we see then that the cou- 
pling between the superconductor and the ferromagnet 
disappears, as does the proximity effect. Furthermore, if 
we "turn off" the spin-orbit coupling in the barrier (i.e. 
setting ciR = 0, so that q + = q~), we see that the spin- 
mixing terms vanish, 5k — SK t — 0, and we recover the 
standard boundary conditions for the proximity effect. 

To describe the superconductor [see |5|] we intro- 
duce a four-component state- vector in the combined spin- 
Nambu space, 
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The Pauli matrices in Nambu space will be denoted by where 
f x , f y and t z , and the unit matrix by tq. 4x4 matrices 
in the spin-Nambu space will be denoted by an inverse 
hat, for instance, 



T p 

f* 



Rp 

R* 



(23) 



where the asterisk denotes complex conjugation. Note 
that 7_ p = 7 p + 7r, therefore, exp(i7_ p ) = — exp(i7 p ), 
as it appears in Tt p and R*_ p - 

Using these notations, we can write the boundary 
conditions for the superconductor in the following com- 



pact form, ma0,#W(- a ) 



'(0)- 



'HO- 



We proceed by first incorporating the right side of this 
boundary condition into the Bogoliubov-de Gennes equa- 
tion for a superconductor, which we symbolically write 
in the formj^ 



2rriB 



5{z- 



T p ^)(0)-i? p ? p *)( 



(24) 



Here, G p (z,z') is the Green function of the Bogoliubov- 
de Gennes equation for a bulk superconductor in the 
half-space, z < —a, which satisfies the von Neumann 

boundary conditions, d z G p (z, z 1 ) = 0. It can be ex- 

z— — a 

pressed in terms of the Green function of an infinite bulk 
s-wave superconductor by the method of mirror images 
and reads, 



G p (z, z') = G p (z - z') + G p {z + z' + 2a), 



where 



P " \F p G*_ p 



(25) 



(26) 



and G p and F p are the normal and Gor'kov Green func- 
tions, respectively. Consequently, the solution to the 
Bogoliubov-de Gennes equation (24) is derived by con- 



voluting the Green function with the boundary term on 
the right side of (24), which is simple due to the delta- 



function in the latter, 



m B 



G p (z + a) Tp^ p F )(0) - R p ^ p s) (~a) 



(27) 

For the purpose of understanding induced supercon- 
ductivity in the ferromagnet we only need to know the 
boundary value of the superconducting wave function at 
z = —a. From (27) we obtain 



1 



m B 1 - {h 2 /m B )g P R P 



g p fjf (f)(0), 



(28) 



9 P = lim G p {z + a)= [ ( ^G(p x ,p y ,p z ). 

z^ — a J ZTT 



(29) 



We note here that ( 28 ) is exact and that we have not used 



any properties of the system at z > up to this point. 
Therefore, the equation above is applicable to any such 
junction with any normal or superconducting material at 
z > (the only constraint is the continuity of derivatives 
at the z — boundary, which may be violated if there 
is a spin Hall effect present. In this case, however, an 
alternative set of boundary conditions can be derived) . 



3. Ferromagnet 



The expression in ( 28 ) , which represents a useful tech- 



nical result of the paper, allows one to formulate a closed 
problem on the (ferromagnetic) normal side. Let us write 
the Schrodinger equation ([3| as 



G(f)- 1 o^(f)(z)=0, 



(30) 



where we have extended the equation into the Nambu 
space. The propagator G p [z, z ) describes a particle of 
the magnetized Fermi liquid in the shell < z < d. Apart 
from the trivial boundary condition ([6| at the hard wall, 
^ together with the solution of Sec. 
give rise to the following constraint, 



IIIB1 and (28) 



Rr> + T„ 



h 2 /ttlb 



7T9p T p 



* P F) (o). 

(31) 



The boundary condition ^ together with (30) and (31) 
form a self-consistent set for z > 0. 

We assume at this point that the spin-orbit coupling 
energy scale is small compared to other relevant energies 
in the problem and keep the spin-orbit parameter cxr fi- 
nite only where we otherwise would get a vanishing effect, 
i.e. in the spin-mixing terms. In all other quantities, we 
set an = 0. This brings the reflection matrix to a simpler 
form proportional to the unit matrix, 



kI — KTQffQ, 



where k is defined in (21 ). In the cxr 



(32) 

limit we have 
e). In the 



2m 



k = q/ tanh (qa), where q - 
limit of a high barrier, k q. 

We now consider the operator denominator in the 
boundary condition (31). If the bulk superconductor 



is deep in the paired state, we may neglect the normal 
Green function and estimate the integrated Green func- 
tion in (29) as g sa ff x (ia y ), with 



A 



A 2 



(33) 



G 



For small £ p = p 2 /(2ms) — /is, it can be estimated as 
/ ss l/(2hvs), where vs is the Fermi velocity in the su- 



perconductor. Therefore, the term in question from (31) 
becomes 



1 



m B 1 



H 2 



9 P 



1 



n 1 + a 2 



[T x {iv y ) - Oil] , (34) 



where a = Kf/rriB- This parameter can be estima ted as 



a ~ \l msUo/mBEp?\ Note that the last term in (34) is 



uninteresting because it does not include any off-diagonal 
contributions in Nambu space. The term slightly renor- 
malizes the boundary conditions for the transverse wave 
function in the ferromagnet. It can therefore be safely 
dropped. 

Incorporating the right side of the boundary condition 



(31) into the Schrodinger equation for the ferromagnet 
(30) we obtain 



-h 2 a8(z) 



2/kto b (1 + a 2 ) 



This is the Bogoliubov-de Gennes equation for the nor- 
mal region with superconducting correlations introduced 
via the boundary term that carries information about all 
Andreev processes. 



function taken proportional to ia y is basis-dependent and 
assumes spin quantization along the z-axis). 

For the ferromagnetic state polarized along z, we ob- 
tain the Bogoliubov-de Gennes equations in the familiar 
form, 



Imp 



- {^F ~ £tr) 



E g {p v + ip x )4>% = 0, (39) 



where the proximity-induced gap or so-called minigap is 



E„ 



—^-KfSKt 



V> (tr) (o) 



(40) 



The boundary value of the transverse wave function can 
be obtained from (37) and reads 



i/; (tr) (0) 



1 



di+e + i/{Kdy 



(41) 



where £ = mpq/ '(msfctr) an d &tr is a solution to the 
eigenvalue problem, tan (k tI d) / (k tl -d) — —ms /(mpqd), 
which determines the spectrum, e tr = h 2 k 2 r /(2mB)- We 
see from Eq. (40) that the proximity-induced minigap 



is extremely sensitive to the actual width of the normal 
region: the smaller the width, the more frequent Andreev 
scattering processes, and the larger the minigapP^ 



C. Proximity-induced superconducting gap 

In order to calculate the superconducting gap induced 
in the ferromagnet, we assume that its magnetization is 
unaffected by the proximity to the superconductor and fix 
it instead of calculating it self-consistently. We also focus 
on a particular sub-band with a given spin polarization 
and ignore intraband scattering (which is justified only 
in the clean case). Furthermore, we assume that the 
ferromagnetic layer or wire is very thin in the z-direction 
and approximate the full solution as 



^W=Xo(m)V (tr) W0(p), 



(36) 



where Xa( m ) is a spinor describing the magnetization di- 
rection that we wish to enforce and tp^ tT \z) is the trans- 
verse envelope wave function that we approximate as a 
solution of the free one-dimensional Schodinger equation, 



d 1 



ZmpEtr 

h 2 



V (tr) (z) = o, 



(37) 



with the unusual boundary conditions, 



m B 
mp 



a z V (tr) (0) = «^ (tr) (0), V (tr) (d) = 0. (38) 



We then introduce solution (36) into (35), multiply it 
by ip^* (z) , average it over the transverse direction, and 
extract the relevant spin component (a word caution here 
is that the standard convention for the Gor'kov Green 



IV. RESISTANCE PEAK DUE TO 
SUPERCONDUCTING FLUCTUATIONS 



We argued in Sec. |II B| that the W electrodes have a 
strong impact on the transport properties through the Co 
wire. We reiterate the experimental fact that when an 
extra W strip is deposited on top of the Co wire [see Fig. 
[ljc)], the normal state resistance of the wire increases 
by nearly 50%. This indicates that when a W electrode 
is deposited onto a nanowire they modify the wire in 
its vicinity and provide a major source of resistance. A 
faithful model of the Co wire may then be such that the 
current i w goes through the two W voltage electrodes. 
If this view is adopted this creates four W-Co interfaces, 
each one similar to the interface considered in Sec. Mil 
With this geometry electron transport through the wire 
is expected to be strongly modified by Andreev physics 
at the W-Co interfaces. 

We now focus on the resistance peak observed near the 
superconducting transition temperature of the W elec- 
trodes. This peak is very reminiscent of similar anoma- 
lous peaks studied in the context of c-axis transport in 
cuprate superconductors^^, magnetoresistance in dirty 
filma^l, as well as magnetoresistance in granular elec- 
tronic systems^. In all these cases, the anomalous peak 
has been explained using the phenomenon of supercon- 
ducting fluctuations ! 41 * 42 ^ This phenomenon is concerned 
with fluctuating Cooper pairs that form while the system 
is in the normal state, either just above the transition 
temperature T c or the critical magnetic field H c . For the 



7 



H = case, the appearance of Cooper pairs above T c 
opens up a new channel for charge transport. Indeed, 
these fluctuating Cooper pairs can be treated as carriers 
of charge 2e with a lifetime given by tql ~ ^/^b(7 1 — T c ). 
This leads to a contribution to the conductivity known as 
the Aslamazov-Larkin (AL) contribution or paraconduc- 
tivity, and gives a positive correction to the Drude con- 
ductivity. There is also an indirect correction known as 
the density of states (DOS) contribution. One of the im- 
portant consequences of these fluctuating Cooper pairs is 
the decrease in single-particle DOS near the Fermi level. 
The idea is that if some electrons are involved in pairing 
they can not simultaneously participate in single-electron 
transport. The DOS contribution therefore gives a neg- 
ative correction to the Drude conductivity. 

As we approach T c of the W electrodes, Cooper pair 
fluctuations grow inside the W electrodes. By virtue of 
the Rashba spin-orbit coupling at the W-Co interfaces 
(due to breaking of mirror symmetry in the longitudinal 
direction at the interface) singlet Cooper pairs inside W 
are converted and triplet fluctuations permeate into the 
Co wire. AL correction to the conductivity arises due to 
the transport of these fluctuating Cooper pairs through 
the wire. However, Cooper pair tunneling across an in- 
terface is strongly suppressed because it is a higher-order 
process in both tunneling and the singlet-triplet conver- 
sion rate, which requires a spin flip process enabled only 
by weak interfacial spin-orbit coupling. This can be seen 



in (17 1 where the off-diagonal elements of the transmis- 
sion matrix are suppressed by the factor Sn t , which is 
proportional to a^. The single electron tunneling and 
the DOS correction, on the other hand, are lower order 
tunneling processes that require no spin flipping. We 
therefore qualitatively expect the DOS correction to be 
parametrically larger than the AL contribution, thus giv- 
ing a negative overall correction to the Drude conductiv- 
ity. This clearly can explain the anomalous upturn in the 
resistance as a function of temperature. 

When the Co wire is replaced by a gold wire, the sit- 
uation is modified. Since gold is not ferromagnetic, s- 
wave pairing correlations are expected to survive over 
a much longer distance inside the wire. In this case, 
Cooper pair tunneling through the W-Au interface is 
much more transparent since it does not require spin flip- 
ping. Therefore, Cooper pairs are more efficient at trans- 
porting charge in this case, and the AL contribution to 
the conductivity is expected to play a much bigger role 
than with a ferromagnetic wire. Here, we purport that 
the absence of the resistance peak in the Au wire can be 
explained if indeed the AL contribution to the conductiv- 
ity is parametrically larger than the DOS contribution, 
thus giving an overall positive correction to the conduc- 
tivity. 

Near the superconducting transition temperature (and 
at H = 0) and in dimensions at or below two, both AL 
and DOS contributions show either algebraic or logarith- 
mic divergent behavior as a function of the distance from 
the transition point. Denoting the divergent part of the 



AL corrections by function fAh(T—T c ) and the divergent 
part of the DOS correction by Jdos(T—T c ) (and ignoring 
the Maki-Thompson contribution), the total supercon- 
ducting fluctuation correction to the normal resistance 
Rq can be schematically written as 



6R{T) 

Rq 



= ^dos/dos(T - T c ) - A Ah f AL (T - T e ), (42) 



where, Aal and Ados & re positive pre-factors. As we 
will show in Sec. |IV A| 



1 



while 



AL 



/DOS 



(T/T c - I) 4 



yfTjT~ c 



(43) 



(44) 



indicating that the AL contribution diverges much more 
strongly than the DOS term as T — >• T+. Therefore, un- 
less the pre-factor Aal is exceedingly small in compar- 
ison to Ados the DOS contribution is always paramet- 
rically smaller than the AL contribution and the con- 
ductivity correction is positive. This leads to a mono- 
tonic decrease in the resistance as the system approaches 
the transition, and no peak will result. However, if we 
have a situation where the AL pre-factor is greatly sup- 
pressed, i.e. Aal "C Ados: the DOS contributions may 
become parametrically larger than the AL counterpart 
for T sufficiently, but not too close to T c + . In this case, 
the resistivity can display an anomalous upturn before 
the AL contribution eventually takes over and the resis- 
tivity makes a precipitous fall as T — > T+. 



A. Microscopic theory 

1. Effective dimensionality of a tungsten electrode for 
fluctuation analysis 

The superconducting fluctuation physics depends crit- 
ically on the dimensionality of the system, with the gen- 
eral trend being that the lower the dimensionality, the 
more pronounced and singular the fluctuation effects arc. 
However, one should exercise care in determining the ef- 
fective dimensionality of a system, as this notion depends 
on a particular effect that is being studied. For example, 
the system may be three-dimensional in terms of single- 
electron diffusion physics, but fall into the category of 
one-dimensional superconductors when it comes to the 
fluctuation analysis. We believe that such is the case 
with the W electrodes that host superconductivity and 
most of the fluctuation physics in the experiment under 



(i) 



(2) 

250nm wide and L y ± ' ~100nm 



study involving L ± 
thick W strips. 

The finiteness of the transverse directions implies that 
whenever we have an integral over a three-dimensional 
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as follows 




FIG. 4: (a) Self-energy diagram describing a correction to the 
single-electron Green function. The wavy line corresponds to 
the superconducting fluctuation propagator, and the shaded 
vertices represent Cooperon vertices, (b) A diagram corre- 
sponding to the lowest-order correction to the transport kernel 
due to superconducting fluctuations in the tungsten electrode. 
As labeled in the figure, top (bottom) solidline represents the 
electron propagator for tungsten (cobalt). 



momentum, its transverse part must be replaced with a 
sum over quantized modes, 



dg± 
2nh 



f(qx) 



1 ^ f ( 27r?m_i_ 



For superconducting fluctuation analysis, the relevant 
function of interest is the fluctuation propagator, which 
appears in the combination 



1 



E 



ii 



'GL 



(45) 



where D = Wpr/3 is the diffusion coefficient and the 
Ginzburg-Landau relaxation time, 



t~gl 



T-T c 



(46) 



tunes the proximity to the transition temperature T c . 
Hereafter, we will focus on a superconducting electrode 
and assume that we are dealing there with a disordered 
superconductor. 



The question of whether or not a particular dimension is 
important reduces to the comparison of the first and last 
term in the square brackets in ( 45 1 . If the former is much 



larger than the latter for any n± ^ the corresponding 
dimension is unimportant and an effective reduction of 
dimensionality occurs. One can define the following char- 
acteristic temperature scale (to be compared with T — T c ) 



Tl(Li 



1 



I 



( ^ v f \ 



(47) 



where i>f is the Fermi velocity, I = v^t is the electron 
mean-free path and we have restored the Boltzmann con- 
stant /cb and the Planck constant H. 

For W, the Fermi velocity can be estimated as v-p ~ 
0.5 • f0 6 m/s and taking the largest of the two transverse 
dimensions Lj_ ~ 250 nm, we find 



Tf(250nm) 



I 



80 K. 



If we assume that the mean-free path in the amorphous 
W strips is of order a nanometer the corresponding tem- 
peratures scale becomes of order T^(250nm) ~ 0.3-1K, 
which is very reasonable and implies that as soon as we 
approach the superconducting transition with (T~T C ) <C 
T c ~ 4.4-5K, we may view the electrode as a one- 
dimensional superconductor. 



2. Fluctuation correction to the density of states 

We now analyze the DOS fluctuation physics on the 
basis of the standard diagrammatic perturbation theory. 
The Cooper-channel correction to the electronic DOS is 
given by 



6u(e) 



1 



/ (2^ Im {£ 2 ( ie ™>P) E ( i£ ™>P)} 



7T J (27r) 



(48) 

where £(ie m ,p) is the self-energy described by the dia- 
gram in Fig. |4][a). It reads 



E(ie m ,p) 



where 



/d d q 

xC 2 (£ m ,fi n -e ra ;q)£(fl„,q), (49) 



Q(ie m ,p) = — 



(50) 



is the Matsubara Green function with £ p = v-p(p — p-p) 
and s m = e m +sgn (e m )/ (2r) and r is the scattering time, 

C(e m ,n n - e m ;q) = - ; - ; ^- ^ (51) 



r Dq 2 + 7a + |2e r , 



is the Cooperon vertex, where we included pair-breaking 
scattering rate 7 S = 2/r s and 9{ ) is the standard Heavi- 
side step function. In the vicinity of the transition point 
the fluctuation propagator reads 



c(n n , q )= 8 ^[Dq 2 

7TI/ 



'GL 



M \ (52) 



where tql = f T ^ T and vq = mpp /2tt 2 is the bare DOS 
for a single spin projection. As per the usual convention, 
n n = 2imT denotes the bosonic Matsubara frequency 
and e m = (2m + l)wT is the fermionic Matsubara fre- 
quency. Finally, note that the physical quantity of in- 
terest, the DOS, must be analytically continued from the 
discrete set of Matsubara frequencies to the continuum of 
real energies, as labelled by the symbol ie m —> e in ( 48 1 . 



where 



Since the Tungsten electrodes are amorphous, we can 
safely assume that they are strongly disordered s-wave 
superconductors, where T c t <C 1 (but of course we also 
assume that we are far from localization, that is E^t 3> 
1). In this case, the three- Green-function block takes the 
especially simple form 

v I d£G(in n - ie m ,q - p)G 2 {is m ,p) 



= -27rw r 2 sgn(e m ), (53) 
where we enforced the constraint e m (e m — Q n ) > in 



(51 1 . This leads to the following expression for the DOS 
d d q ( d 

where S R (q,e) — S(q,ie m —> e) represents the 
analytically-continued (retarded) Matsubara sum 



Me) = 



Re^(q,e), (54) 



5(q,e m > 0) 



(55) 



E 

n=— oo 



Dq 2 



2ttT c 



Dq 2 +ls 
2ttT c 



1 2m + 1 - n\ 



This sum can be calculated exactly in terms of the 
digamma function, 



oo 

7 + E [{n + l)' 1 - {n + zY 
n=0 



(56) 



which is analytic everywhere except z = 0, — 1, —2, .... It 
is convenient to separate the sum ( 55 ) into two pieces, 

Y^ n =-oo ' ' ' an d J2™=o "'• For the nrs t term analytic 
continuation reduces to replacing e m — ¥ —ie. For the 
second term it becomes possible after noticing the "re- 
flection property" , where m—n can be replaced by n— m. 
This leaves the sum from n — to n = m unchanged but 
the denominator "positively-defined" and ready for ana- 
lytic continuation. Finally, the asymptotic form of DOS 
in the limit of {e, Tq^, 7 s } «T ~ T c becomes 



S*(q,e) 



(27rr c ; 



(Dq 2 



' GL, 



(2Dq 2 + 7s + T G l 



2ie) 



(57) 



Let us now focus specifically on the correction to the 
DOS of electronic states with spin-up at the Fermi level, 
i.e. e = 0. The remaining elementary integrals in (54) 
can now be easily calculated and we find 



6v t (0,T) 



TV- V2 



1 + 2a{T) 



P 2 A 128 (T-T c )5/2 Q 3( T) [ 1 + a ( T) ] 



a(T) 



4(T-T c )t s \ 



(59) 



Here, A is the effective area of the fluctuating super- 
conductor in the region where the fluctuations are stud- 
ied (we expect it to be related to the dimensions of the 
superconductor-nanowire contact and, consequently, ex- 
pect A to be smaller than the cross-sectional area of the 
wire). Notice that in the absence of pair breaking and/or 
far from the transition, where (T — T c )t s S> 1, |59| ) re- 
duces to a constant a = l/y2 and the DOS at the Fermi 
level acquires a very sharp temperature dependence as 
follows: 

**t(0,T) cx -(T-T c )- 5 / 2 , if (T-T c )r s »l. (60) 

In the opposite regime of strong pair-breaking scattering 
or in the immediate vicinity of the transition, we find 

<J*t(0,T) cx -(T-Te)- 1 / 2 , if (T-T c )r s <l. (61) 



3. Fluctuation correction to the contact resistance 

The result for the DOS at the Fermi level provides a 
useful insight into the physics near the transition and 
illustrates that, as a precursor to the global pairing tran- 
sition, electrons are actively swept from the vicinity of 
the Fermi level due to the formation of the fluctuat- 
ing Cooper pairs. However, the stand-alone quantity 
Sv(0, T) is not extremely useful for comparison with ex- 
periment, as it is actually not directly measurable. What 
is measured in experiment is resistance, which is an in- 
tegral quantity that includes excitations with different 
energies and that has contributions from both single- 
electron transport across the W-Co junction and pair 
transport. As found in Sec. |III[ the latter is suppressed 
strongly far from the transition due to the very small pair 
tunneling probability, which requires spin-orbit-assisted 
spin flips. Single-electron transport on the other hand 
does not rely on any spin-orbit coupling, and spin-up 
electrons can tunnel freely from the electrodes into the 
ferromagnet. Hence, there is a regime close to the tran- 
sition (but still far enough from the immediate vicin- 
ity, where the Cooper pairs eventually take over) where 
the single-electron tunneling dominates and is already 
strongly suppressed by fluctuations. This results in the 
upturn in the resistance as T — > T r + . 

If we now focus entirely on single-electron transport 
and disregard the pair-breaking scattering, one would 
tend to conclude (incorrectly) that the tunneling resis- 
tance acquires a contribution proportional to that in ( 60 1 . 
This, however, is not so and the correction to the con- 



(58) 



ductance is much weaker. This is because the tunneling 
electrons involve not only those precisely at the Fermi 
level, but all electron excitations within the shell of en- 
ergies E ~ (Ep — T,Ep + T). Hence, what matters is 
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a redistribution of the DOS beyond the shell of energies 
that participate in transport (any redistribution within 
the shell is essentially not observable). 

All these phenomena are automatically accounted for 
within the standard diagrammatic theory of transport. 
Here, we proceed nearly in a one-to-one correspondence 
with Ref. 1561 The tunneling current is derived from the 
transport kernel Q(u>i) shown in Fig. |4^b), 



I(V) = -eImQ R (u) 



j=eV 



(62) 



where Q R (u) is analytically-continued transport kernel 
and V is the voltage across the tunneling contact. In 
the setup under consideration the Matsubara transport 
kernel reads explicitly 

Q{ui) = -T 2 J2 l^k,k'| 2 ^(em,k)g Co (e m + W/ ,k') 

£ m ;0„;k,k' 

^— j£(O n ,q)C 2 (e m ,0„-£ m ;q) (63) 

x g w (Q n - e m ,q- k), 

where ik,k' is the tunneling amplitude between two mo- 
mentum states, Gw/Co(Sm, k) is the electron Green func- 
tion for the W electrode/Co nanowire and all other quan- 
tities have been defined in Sec. IV A 2 Following Ref. 
l56l we introduce a normal state resistance of the W-Co 




FIG. 5: Correction to the differential conductance from 
superconducting fluctuations as a function of dimensionless 
voltage eV/2irT (see (65 1). The three curves correspond to 
different temperatures: t := (T — T c )r = 0.01 for the solid 
line; t = 0.05 for the dashed line; and t = 0.1 for the dotted 
line. 



where Q(z) is the Riemann zeta-function and £(3) ~ 
1.202. The result in (66) corresponds to a sharp upturn 



in the resistance upon approaching the superconducting 
transition of the electrode. 



junction R 



(o) 

W/Co 



and obtain 



7rT 2 

Q( u >l) = 7m — sgn(e m ) sgn (e m +w ; ) (64) 

Ze n W/Co e m ,Q n 



x 6 [e m (e m - fi n )] 



dq 



£(f2„,q) 



2nA (Dq2 + \2e m -n n \y 



Repeating the dimensionality-independent summations 
as in Ref. [SHI (note that the clean case^ becomes tech- 
nically equivalent to the dirty case if we notice that the 
frequency-dependence of the three-Green-function block 
in the reference is identical to that in the two Cooper- 
ons that appear in ( |64[ )) and evaluating the remaining 
integral over momentum q we obtain the non-linear I-V 
dependence 



I(V) 



T 
8^ 



1 



R 



(0) 

W/Co 



Tm-0' 



ieV 



2ttT 



V(T-T c )t 



(65) 



where V'(') is the logarithmic derivative of the gamma- 
function. Now focusing on the linear response regime, 
^\v=o, we find the leading fluctuation correction to the 
W-Co contact resistance as follows, 



SR 



(DOS) 
W/Co 



7C(3) [to 1 



(o) 

W/Co 



R 



4tt V 2 pjU y/(T-T c ) T ' 



(66) 



4- Estimate of tunneling resistance due to fluctuating pairs 

As emphasized throughout this section, the fluctuat- 
ing Cooper pairs that appear in the electrodes do not 
aid with transport through the ferromagnetic nanowire 
initially because of their poor ability to tunnel into it. 
However, the experimental fact that the wire as long as 
a micrometer does become superconducting ensures that 
the pairs eventually are able to tunnel. Sec. [Hi] provides a 
microscopic picture of the Andreev reflection/boundary 
physics that presumably makes this possible. Therefore, 
we expect that as the temperature is tuned down to the 
closest vicinity of the transition, the Cooper pair tunnel- 
ing takes over the effect of the suppression of the DOS 
and the upturn in the resistance crosses over to the down- 
turn going through a peak, as observed in experiment. 

We notice here in passing that the height of the peak 
may provide a valuable insight into the competition of 
the two phenomena and may potentially become a means 
to measure the boundary spin-orbit coupling that is cru- 
cial for the proximity-induced p-wave superconductivity 
in the wire, the phenomenon of major interest. How- 
ever, we leave the complicated microscopic theory of spin- 
orbit-assisted fluctuating pair tunneling for future studies 
and in this section only extract the leading temperature 
dependence of this AL type correction. It can be done 
on the basis of the Ambegaokar-Baratoff formula that 
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yie 



■IdiP 



di V/Co 



-5(o!ij) / de 



cosh 2 (^) 



(67) 



wher e <5t/(e) is the suppression of the DOS obtained in 
Sec. IV A 2 and we kept a small coefficient 6 (an) that 
includes spin-orbit suppression of the pair-tunneling as 
obtained in Sec. Mil 

Note that the integral of the DOS over all energies van- 
ishes identically /_ de8v(e) = 0, (the physical interpre- 
tation being that the total electron density is conserved 
and the electrons can only be redistributed across differ- 
ent energies) and a non-zero correction to the resistance 
in the first order appears only due to the modulating 
function, cosh -2 (h^), and consequently is much weaker 



(see (66l) than that in the DOS (see ([60])), For the pair 
tunneling, the situation is different as de8v 2 (e) is 
not only non-zero, but is very singular function near the 
transition. This singularity determines the scaling of the 
AL correction (67) as 



dR W/Co K " 



S(or) 
(T - T c ) 4 



(68) 



V. TOPOLOGICAL SUPERCONDUCTIVITY 
AND MAJORANA FERMIONS 

A. Relation to the experiment in Ref. 1211 

The canonical Kitaev modeP^, which is the simplest 
prototype of a one-dimensional topological supercon- 
ductor, involves one fermion species hopping in a one- 
dimensional chain and subject to a prescribed p-wave 
pairing field on the nearest-neighbor bonds. The end 
points of the chain host single Majorana excitations that 
are of key interest. One can consider Af replicas of 
the Kitaev-Majorana model with different p-wave pairing 
fields for each species but with no mixing between them, 
i.e. no intcrband pairing. If Af is even, the end Ma- 
joranas are unstable against various perturbations and 
generally hybridize into relatively uninteresting finite- 
energy boundary states. In contrast, if N is odd, the 
system can then host one Majorana zero-energy state at 
each end. Recently, theoretical works have addressed a 
multi-channel generalization of Majorana end states in 
quasi-one-dimensional structures with Rashba spin-orbit 
coupling! 57 " 60 ! The works show that the Majorana end 
states are realized in some parameter regime as long as 
an odd number of transverse subbands are occupied. 

The well-known difficulty in realizing the regime where 
(# of end Majoranas) mod 2 = 1 is of course re- 
lated to the fact that the electrons have spin and 
consequently come in pairs. Therefore, according to 
Kitaev^, one necessarily, but not sufficiently, needs to 
break time-reversal symmetry in the one-dimensional 



superconductor to render it topological. The notable 
and fascinating proposals, where this regime was the- 
oretically shown to be possible, include a heterostruc- 
ture involving 5- wave s uperconductor/ spin-orbit-coupled 
wire/ferromagneP°J22l, a similar heteros tructu re with- 
out a ferromagnet but in a magnetic field 31 37 , topolog- 
ical insulator/superconductor in a fiel d 28 * 40 *, and half- 
metal/superconductor with spin-orbit coupling in the 
latteM 

If our interpretation of the experimental data in Ref. 
l2"Tlis correct, we see that this work potentially possesses 
all necessary ingredients for the realization of the Majo- 
rana end states; it has a ferromagnetic crystalline wire 
in which p-wave superconductivity has been induced. 
The A^-replica Kitaev-Majorana model discussed here 
may apply to a thin narrow Co film deposited on top 
of a three-dimensional W superconductor. The Co film 
should be highly confined in the direction normal to the 
Co-W interface so that only one channel is occupied 
in that direction, and the width of the strip should be 
smaller than the p-wave coherence length. Since Co is not 
a half-metal the question is whether a ferromagnetic wire, 
with both majority and minority carriers, can still host 
an odd number of end Majorana fermions. This, in prin- 
ciple, may be possible, as what we need is not necessarily 
to eliminate completely a spin component, but merely 
to make the two components different from each other 
such that the total number of occupied subbands is odd, 
i.e. (A/f + A/4.) mod 2 = 1. In the simplest model with 
no inter-subband mixing (considered here), this would 
imply topological superconductivity. However, we men- 
tion that a single-species p-wave superconductor, more 
akin to the original proposal by Kitaev 2 - 7 -', may be possi- 
ble if the thin film Co wire is replaced by a similar film 
made of a half-metal such as Cr02. Notably, experiments 
have recently shown that Cr02 can accommodate long- 
ranged p-wave superconductivity when it is proximity- 
coupled to a conventional superconductor.^ Majorana 
bound states should then be realized at the ends of such 
a half-metal wire when odd number of transverse sub- 
bands are occupied.^3 



B. Realizing Majorana end states using a 
ferromagnetic semiconductor / ferromagnetic 
superconductor heterostructures 

Ref. [21] and the discussion above suggest an even sim- 
pler experimental setup for topological superconductiv- 
ity. Indeed, realizing a one-dimensional topological su- 
perconductor using spinful fermions requires lifting the 
double degeneracy imposed by time-reversal symmetry. 
In principle, this can be achieved by proximity-inducing 
a Z eema n gap^l or by applying an external magnetic 
fiel d 31 ' 36 '. An alternative approach to realizing a one- 
dimensional topological superconductor is to deposit a 
ferromagnetic semiconductor wire on top of a ferromag- 
netic superconductor^^. A particularly attractive can- 
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didate ferromagnetic semiconductor is europium oxide 
(EuO) , whi ch is known to possess nearly spin-polarized 
bands EuO becomes ferromagnetic below 70K under 
ambient pressure and the Curie temperature is known 
to increase with pressure reaching 200K under 1.5 x 10 5 
atmospheres Since its integration with with Si and 
GaN^l, EuO has garnered much attention for its po- 
tential use in spintronic applications. Ferromagnetic su- 
perconductors are materials that exhibit intrinsic coexis- 
tence of superconductivity and ferromagnetism where the 
same electrons are believed to be superconducting and 
ferromagnetic simultaneously. Perhaps the most well- 
known experimental realizations of ferromagnetic super- 
conductors are the Uranium-based compounds. Follow- 
ing the first experimental realization a decade ago there 
are now four such compoundiP^, two of which exhibit 
the phenomena at ambient pressures^'. We propose that 
depositing URhGe electrodes on the EuO wire, in the ar- 
rangement shown in Fig. [lja), would be conceptually 
the simplest structure to realize Majorana fermions that 
would require neither topological insulators nor control 
over spin-orbit couplings. Application of a magnetic field 
can further stabilize both the ferromagnetism and super- 
conducting phase. 

VI. SUMMARY 

In this work, we study temperature-dependent trans- 
port properties of a ferromagnetic cobalt nanowire 
proximity-coupled to superconducting electrodes. The 
work is largely motivated by a recent experiment,^ in 
which long-ranged superconducting proximity effect was 
observed in single-crystal Co nanowires coupled to con- 
ventional superconductors. We focus particular atten- 
tion to wires where the transition to superconductiv- 
ity is preempted by a large and sharp resistance peak 
near the transition temperature of the electrodes. We 
explore the possibility of inducing spatially-odd p-wave 
correlations in the wire via Rashba spin-orbit coupling 
in the wire-superconductor interface. We then theoret- 
ically study the physics of an s-wave-superconductor- 
ferromagnet tunnel junction in the presence of an in- 
terfacial Rashba spin-orbit coupling, and derive an ex- 
pression for the induced p-wave minigap in terms of the 
microscopic parameters of the contact. 

In the second half of the work, we show how the anoma- 
lous resistance peak can be explained within the theory 
of superconducting fluctuation corrections to conductiv- 



ity. In particular, we develop a microscopic theory for the 
superconducting fluctuation corrections in the tungsten 
electrodes and discuss how these corrections may lead to 
corrections in the measured resistance. In the last part 
of the work, we discuss in detail the possible relevance of 
the experiment to topological superconductivity and pro- 
pose a related hybrid system, which provides the simplest 
physical realization of topological superconductivity and 
localized Majorana modes. 

We note in conclusion that the spin-orbit-assisted 
proximity-effect scenario considered here is different from 
the diffusive ferromagnet case, where the likely mech- 
anism for proxi mity e ffect is believed to be the odd- 
frequency pairing|22122l The main argument for the odd- 
frequency superconductivity is that any anisotropic pair- 
ing would decay fast into a disordered metal on a length- 
scale of order a mean- free-path, while the odd-frequency 
isotropic pairing is immune from non-magnetic static dis- 
order. This result follows from the Usadel equation if we 
assume that the non-s-wave part of the disorder-averaged 
condensate wave-function is small. Then upon lineariza- 
tion of the Usadel equation, further exponential decay of 
the disorder-averaged p-wave condensate wave-function 
inevitably follows. We note however that in certain ge- 
ometries (e.g., if the mean- free path is much larger than 
the cross-section of the wire, but much smaller than its 
length), there is no reason to assume that the anisotropic 
component is small locally near the contacts and conse- 
quently the conventional derivation of the Usadel equa- 
tions from the Eilenberger equations breaks down in 
these regions. With such assumption about the bound- 
ary conditions, the correlator of p- wave condensate wave- 
functions is not expected to decay exponentially at T = 0, 
in sharp contrast to the exponential decay o f the dis order- 
averaged p-wave condensate wave-functionPSEJESI These 
mesoscopic fluctuation effects will be considered in detail 
elsewhere PS 



Acknowledgments 

V. G. would like to thank K. Efetov, S. Das Sarma, Y. 
Oreg, J. Paglione, G. Tkachov and A. Varlamov for dis- 
cussions. S. T. would like to thank J. Wang for sharing 
their experimental data, and J. R. Anderson for infor- 
mative discussions on the cobalt band structure. S. T. 
would also like thank J. D. Sau for stimulating discus- 
sions at an early stage of the work. S. T. was supported 
by DOE and V. G. by DOE-BES (DESC0001911). 



1 S. S. Saxena, P. Agarwal, K. Ahilan, F. M. Grosche, R. 
K. W. Haselwimmer, M. J. Steiner, E. Pugh, I. R. Walker, 
S. R. Julian, P. Monthoux, G. G. Lonzarich, A. Huxley, I. 
Sheikin, D. Braithwaite and J. Flouquet, Nature 406, 587 
(2000). 

2 D. Aoki, A. Huxley, E. Ressouche, D. Braithwaite, J. Flou- 



quet, J. -P. Brison, E. Lhotel and C. Paulsen, Nature 413, 
613 (2001). 

3 N. T. Huy, A. Gasparini, D. E. de Nijs, Y. Huang, J. C. P. 
Klaasse, T. Gortenmulder, A. de Visser, A. Hamann, T. 
Gorlach and H. v. Lohneysen, Phys. Rev. Lett. 99 67006 
(2007). 



13 



4 T. Akazawa, H. Hidaka, T. Fujiwara, T. C. Kobayashi, E. 
Yamamoto, Y. Haga, R. Settai and Y. Onuki, J. Phys.: 
Condens. Matter 16, L29 (2004). 

5 C. Pfleiderer, M. Uhlarz, S. M. Hayden, R. Vollmer, H. v. 
Lohneysen, N. R. Bernhoeft and G. G. Lonzarich, Nature 
412, 58 (2001). 

6 L. Bauernfeind, W. Widder and H. Braun, Physica C 254, 
151 (1995). 

7 C. Bernhard, J. L. Tallon, C. Niedermayer, T. Blasius, 

A. Golnik, E. Briicher, R. K. Kremer, D. R. Noakes, C. 
E. Stronach and E. J. Ansaldo, Phys. Rev. B 59, 14099 
(1999). 

8 D. J. Pringle, J. L. Tallon, B. G. Walker and H. J. Trodahl, 
Phys. Rev. B 59, 11679 (1999). 

9 V. G. Hadjiev, A. Fainstein, P. Etchegoin, H. J. Trodahl, 
C. Bernhard, M. Cardona and J. L. Tallon, Phys. Stat. Sol. 
(b) 211, R5 (1999). 

10 S.K.Sinha, G.W.Crabtree, D.G.Hinksand, and H.Mook, 
Phys. Rev. Lett. 48, 950 (1982). 

11 J.W.Lynn, G.Shirane, W. Thomlinson, R.N.Sheltonand, 
and D.E.Moncton, Phys. Rev. B 24, 3817 (1981). 

12 J.W. Lynn, J.A. Gotaas, R.W. Erwin, R.A. Ferrell, J.K. 
Bhattacharjee, R.N. Shelton, and P. Klavins, Phys. Rev. 
Lett. 52, 133 (1984). 

13 C.-K. Hsu, D. Hsu, C.-M. Wu, C.-Y. Li, C.-H. Hung, C.-H. 
Lee and W.-H. Li, J. Appl. Phys. 109, 07B528 (2011). 

14 D. A. Dikin, M. Mehta, C. W. Bark, C. M. Folkman, C. 

B. Eom and V. Chandrasekhar, arXiv:1103.4006vl. 

15 J. A. Bert, B. Kalisky, C. Bell, M. Kim, Y. Hikita, H. Y. 
Hwang and K. A. Moler, Nature Physics 7, 767 (2011). 

16 A. I. Buzdin, Rev. Mod. Phys. 77, 935 (2005). 

17 M. Giroud, H. Courtois, K. Hasselbach, D. Mailly and B. 
Pannetier, Phys. Rev. B 58, R11872 (1998). 

18 V. Petrashov, I. Sosnin, I. Cox, A. Parsons and C. Troadec, 
Phys. Rev. Lett. 83, 3281 (1999). 

19 V. Pena, Z. Sefrioui, D. Arias, C. Leon, J. Santamaria, M. 
Varela, S. J. Pennycook and J. L. Martinez, Phys. Rev. B 
69, 224502 (2004). 

20 R. S. Keizer, S. T. B. Goennenwein, T. M. Klapwijk, G. 
Miao, G. Xiao and A. Gupta, Nature 439, 825 (2006). 

21 J. Wang, M. Singh, M. Tian, N. Kumar, B. Liu, C. Shi, J. 
K. Jain, N. Samarth, T. E. Mallouk and M. H. W. Chan, 
Nature Physics 6, 389 (2010). 

22 F. S. Bergeret, A. F. Volkov and K. B. Efetov, Rev. Mod. 
Phys. 77, 1321 (2005). 

23 F. S. Bergeret, A. F. Volkov and K. B. Efetov, Phys. Rev. 
Lett. 86, 4096 (2001). 

24 V. L. Berezinskii, JETP Lett. 20, 287 (1974). 

25 M. Eschrig, J. Kopu, J. C. Cuevas and G. Schon, Phys. 
Rev. Lett. 90, 137003 (2003). 

26 M. Eschrig and T. Lofwander, Nature Physics 4, 138 
(2008). 

27 A. Yu. Kitaev, Phys.-Usp. 44, 131 (2001). 

28 L. Fu and C. L. Kane, Phys. Rev. Lett. 100, 096407 (2008). 

29 Y. Tanaka, T. Yokoyama, and N. Nagaosa, Phys. Rev. 
Lett. 103, 107002 (2009). 

30 J. D. Sau, R. M. Lutchhyn, S. Tewari and S. Das Sarma, 
Phys. Rev. Lett. 104, 040502 (2010). 

31 R. M. Lutchyn, J. D. Sau and S. Das Sarma, Phys. Rev. 
Lett. 105, 077001 (2010). 

32 Y. Oreg, G. Refael and F. v. Oppen, Phys. Rev. Lett. 105, 
177002 (2010). 

33 V. Shivamoggi, G. Refael and J. E. Moore, Phys. Rev. B 
82, 041405(R) (2010). 



34 P. A. Ioselevich and M. V. Feigel'man, Phys. Rev. Lett. 
106, 077003 (2011). 

35 S. B. Chung, H.-J. Zhang, X.-L. Qi, S.-C. Zhang, Phys. 
Rev. B 84, 060510 (2011). 

36 P. A. Lee, arXiv:0907.2681. 

37 J. Alicea, Phys. Rev. B 81, 125318 (2010). 

38 M. Duckheim and P. Brouwer, Phys. Rev. B 83, 054513 
(2011). 

39 S. Tewari, J. D. Sau and S. Das Sarma, Ann. Phys. (NY) 
325, 219 (2010). 

40 A. Cook and M. Franz, Phys. Rev. B 84, 201105(R) (2011). 

41 A. I. Larkin and A. A.Varlamov, In Superconductivity: 
Conventional and Unconventional Superconductors, ed. K. 
H. Bcnnemann and J. B. Ketterson, (Springer, Berlin, 
2004), Ch. 10, p. 369. 

42 A. A.Varlamov, G. Balestrino, E. Milani and D. V. Li- 
vanov, Adv. Phys. 48, 655 (1999). 

43 J. R. Anderson and J. E. Schirber, J. Appl. Phys. 52, 1630 
(1981). 

44 F. Batallan, I. Rosenman and C. B. Sommers, Phys. Rev. 
B 11, 545 (1975). 

45 C. M. Singal and T. P. Das, Phys. Rev. B 16, 5068 (1977). 

46 G. J. McMullan, D. D. Pilgram and A. Marshall, Phys. 
Rev. B 46, 3789 (1992). 

47 J. W. A. Robinson, S. Piano, G. Burnell, C. Bell and M. 
G. Blamire, Phys. Rev. Lett. 97, 177003 (2006). 

48 F. Zhou and B. Spivak, cond-mat/9906177v3. 

49 A. F. Volkov, P. H. C. Magnee, B. J. van Wees and T. M. 
Klapwijk, Physica C 242, 261 (1995). 

50 G. Tkachov, Physica C 417, 127 (2005). 

51 L. B. Ioffe, A. I. Larkin, A. A. Varlamov and L.Yu, Phys. 
Rev. B 47, 8936 (1993). 

52 G. Balestrino, M. Marinelli, E. Milani, A. A. Varlamov and 
L. Yu, Phys. Rev. B 47, 6037 (1993). 

53 G. Balestrino, E. Milani and A. A. Varlamov, Physica C 
210, 386 (1993). 

54 V. M. Galitski and A. I. Larkin, Phys. Rev. B 63, 174506 
(2001). 

55 I. S. Beloborodov, K. B. Efetov, A. V. Lopatin, V. M. 
Vinokur, Rev. Mod. Phys. 79, 469 (2007). 

56 C. Di Castro, R. Raimondi, C. Castellani, and A. A. Var- 
lamov, Phys. Rev. B 42, 10211 (1990). 

57 A. C. Potter and P. A. Lee, Phys. Rev. Lett. 105, 227003 

(2010) . 

58 A. C. Potter and P. A. Lee, Phys. Rev. B 83, 094525 

(2011) . 

59 R. M. Lutchyn, T. D. Stanescu and S. Das Sarma, Phys. 
Rev. Lett. 106, 127001 (2011). 

60 T. D. Stanescu, R. M. Lutchyn and S. Das Sarma, Phys. 
Rev. B 84, 144522 (2011). 

61 P. G. Steeneken, L. H. Tjeng, I. Elfimov, G. A. Sawatzky, 
G. Ghiringhelli, N. B. Brookes and D.-J. Huang, Phys. 
Rev. Lett. 88, 047201 (2002). 

62 A. Schmehl, V. Vaithyanathan, A. Herrnberger, S. Thiel, 
C. Richter, M. Libcrati, T. Heeg, M. Rockerath, L. F. 
Kourkoutis, S. Miihlbauer, P. Boni, D. A. Muller, Y. 
Barash, J. Schubert, Y. Idzerda, J. Mannhart and D. G. 
Schlom, Nature Mater. 6, 882 (2007). 

63 N. M. Souza-Neto, D. Haskel, Y.-C. Tseng and G. Laper- 
tot, Phys. Rev. Lett. 102, 057206 (2009). 

64 F. Zhou and B. Spivak, Phys. Rev. Lett. 80, 5647 (1998). 

65 V. M. Galitski and A. I. Larkin, Phys. Rev. Lett. 87, 
087001 (2001). 

66 S. Takei and V. Galitski, to be published. 



14 



